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After recalling briefly the connection between spontaneous symmetry breaking 
and off-diagonal long range order for models of magnets a general proof of sponta- 
neous breaking of gauge symmetry as a consequence of Bose-Einstein Condensation 
is presented. The proof is based on a rigorous validation of Bogoliubov's c-number 
substitution for the k = mode operator a . 
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1. Introduction 

The connection between Bose-Einstein Condensation (BEC) and spontaneous gauge 
symmetry breaking (GSB) is often taken for granted, usually with the following argument 
that appeals to Bogoliubov's pioneering work pQ: BEC means that in the thermodynamic 
limit the density V^ 1 (aQa ) of particles in the mode of lowest momentum k = is 
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different from zero. The spontaneous breaking of gauge symmetry, on the other hand, 
is associated with the non-vanishing of V _1 |(a )| 2 in the thermodynamic limit. Here 
it is important that the expectation value is understood in the sense of Bogoliubov's 
'quasi averages', i.e., with a symmetry breaking field that is taken to zero after the 
thermodynamic limit. By the Cauchy-Schwarz inequality, we always have |(a )| 2 < 
( a o a o)- Thus GSB implies BEC. On the other hand, if BEC holds, the macroscopic 
occupation of the zero mode is claimed to justify the replacement of the operators a 
and <2g by c-numbers and hence equating |(a )| 2 with (a,Qa ). Therefore BEC and GSB 
are equivalent. 

A rigorous justification of of Bogoliubov's c-number substitution is not a trivial mat- 
ter, however, and the reasoning above is clearly not a mathematical proof. In a classic 
paper [2], that is less known in the physics community than it deserves, Ginibre es- 
tablished the validity of the substitution as far as the computation of the pressure in 
the thermodynamic limit is concerned. As Ginibre pointed out, it is not necessary to 
assume macroscopic occupancy of the zero mode (although the substitution is only a 
useful computational tool if this is the case.) In a recent paper [3j (see also [4], Appendix 
D) we presented a simple and general proof of the exactness of c-number substitutions 
in bosonic Hamiltonians in the thermodynamic limit, using the Berezin-Lieb inequality 
[3 El [7] (which was only derived after Ginibre's paper appeared). We then proved the 
equivalence of GSB and BEC, where the expectation values are defined as quasi aver- 
ages in the sense of Bogoliubov, and moreover that BEC in the standard sense without 
a symmetry breaking field implies BEC in the quasi average sense and hence GSB. This 
last statement, although intuitively plausible [8], does not seem to have been rigorously 
verified before. In brief: 

BEC -> (BEC) qa ^ GSB. (1) 

The present contribution to the Proceedings of the 21st Max Born Symposium gives an 
account of the proof. For comparison and for pedagogical reasons we start by recalling 
briefly the connection between spontaneous symmetry breaking and off-diagonal long 
range order for a model of a ferromagnet. An important difference between this case 
and that of a Bose gas is that the symmetry breaking term in the latter case does not 
commute with the Hamiltonian. 

Another proof of (BEC) qa <-* GSB, based on Ginibre's original results about the c- 
number substitutions, was given by Siito in [9]. A rather different approach, with focus 
on operator algebraic aspects and correlation inequalities, can be found in a 1982 paper 
by Fannes, Pule and Verbeure [10], see also a recent generalization to nonhomogeneous 
systems by Pule, Verbeure and Zagrebnov [TT] . In these papers the equilibrium states on 
the Weyl algebra for the infinite system and their decomposition into symmetry breaking 
states are studied. Some technical assumptions about these states are required in [lOj 
and [TT] but in Remark 2.1 in [TT] it is claimed that the conclusions are valid under 
weaker assumptions. The algebraic approach sheds a different light on the problem and 
is complementary to our method where essentially the only prerequisites are the existence 
of the thermodynamic limit of the grand canonical pressure and average density. 

Finally we remark that the Berezin-Lieb inequality used in our approach gives explicit 
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error bounds that allow us to go beyond |2J and make c-number substitutions for many 
k-modes at once, provided the number of modes is lower order than the particle number. 
This method, that was used in a different context in [12] . has recently proved to be a 
powerful tool for solving other problems about the Bose gas [T31 [H] . 

2. Symmetry Breaking and Off-diagonal Long Range Order in Magnets 

The concept of spontaneous symmetry breaking means quite generally that a sym- 
metry of the Hamiltonian or Lagrangian of a system is not present in the state under 
consideration (usually a ground state or a thermal equilibrium state). A Heisenberg 
ferromagnet with a Hamiltonian of the form 

H = Hq — BM with H = - J xyS x ■ S y and M = ^ S ( x 3) (2) 

provides a well known example to illustrate this phenomenon. Here A C Z d is a lattice 
cube of finite size | A| and periodic boundary conditions are assumed. The spin operators 
S x = (SP,S^,Si 3) ) with [SW.SfW] = iS^ etc act on a copy of C 2s+1 (with s e 
{2,1,...}) at each lattice site and operators on different sites commute. The coefficients 
J x y > measure the interactions between different sites and B denotes the strength 
of a magnetic field in the 3-direction. If B = the Hamiltonian is invariant under the 
transformations S x — > RS, R £ 0(3). The last term in the Hamiltonian explicitly breaks 
this symmetry if B =/= 0. 

The thermal equilibrium state at temperature T in the finite lattice is defined by 

TrAe~ H / kT 

{A)b,t,K = Tve _ H/kT (3) 
and the magnetization, i.e., the thermal average of the magnetic moment per spin, is 

m A (B,T) = \A\- 1 (M) BiT>A . (4) 
The Gibbs free energy per spin, 

g A (B,T) = -lAl^fcTlnTre-^ 7, (5) 
is a concave function of B with 

m A (B,T) = — . (6) 

Since T may be regarded as fixed in the following we shall drop it in the notation, and 
in order to avoid irrelevant factors, choose energy units so that kT = 1. 
For a finite A, m A {B) is continuous in B, and for B = we always have 



m A (0) = 



(7) 
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(3) 

because the state is invariant under a unitarily implemented rotation that takes S x — > 

(31 

— Sx . On the other hand, the thermodynamic limit of the magnetization 

m(B) = lim m A (B) (8) 

A— Z<* 

may break the symmetry for sufficiently low T, i.e., it is possible that 

m± = lim m(B) ^ 0. (9) 

Since the Hamiltonian is invariant under simultaneous flipping of the spins and B we 
always have m_ = — m + , so if m± 7^ the Gibbs free energy in the thermodynamic 
limit, g(B) = lim^^d g\(B), is not differentiable at B — 0. As a concave function it is, 
however, differentiable almost everywhere. 

We may also consider the average of the square of the magnetic moment per spin, 

i.e., 

ml(B) ((M/\A\f) BTA = (|A|- 2 J] S 3 x S 3 y ) BTA (10) 

x,y£ A 

with a corresponding thermodynamic limit m 2 (B). By the Cauchy-Schwarz inequality 
we always have 

m 2 (B) > (m(B)) 2 . (11) 

Hence, spontaneous symmetry breaking in the sense that lim B ^ + m{B) =/= 0, implies 
off- diagonal long range order in the sense that 

lim m 2 (B) > 0. (12) 

It is important to note, however, that it is a priori not clear that this implies m 2 (0) > 
for reasons explained further below. 

On the other hand, a general argument due to Griffiths [15] shows that 

to 2 (0) < lim m 2 (B)= lim (m(B)) 2 . (13) 

The argument is based on two facts. The first is the commutativity of the symmetry 
breaking term, BM, with the other part, Hq, of the Hamiltonian. This allows a simul- 
taneous diagonalization of the two parts. The second fact is Griffith's lemma that we 
recall for completeness (for a proof see [15] or [16]): 

Griffith's Lemma[15j Let \x n be a sequence of probability measures on M. such that 

f(y) = lim -In / e xy dfi n {x) (14) 

exists for all y in an interval around 0. Define 

a ± = lim MlM, (15 ) 
v^o+ y 
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Then there is a c < 1 such that, for any e > 0, 

/ dtl n (x) = 1 - 0(c n ). (16) 

J (a- —e)n 



In other words: The scaled measure, d/i„(n£), is supported in a_ < £ < a+, with 
exponential accuracy, as n — > oo. 

To apply this lemma to the situation under consideration we take x = B, f(B) = 
g(Q) — g(B), n — |A| and d// n (ra£) the probability distribution at B = for the magnetic 
moment per spin, £ = m = |A| _1 M. Then a± = m±. Moreover, / is a convex function of 
B so the derivative with respect to B is monotonously increasing. We can now conclude 
two things: 1) At points where / is differentiable, i.e., for almost all B ^ 0, the probability 
distribution of the magnetization is a delta function concentrated at m(B) and hence 
m 2 (B) = (m(-B)) . 2) The probability distribution for the magnetization at B = in 
the thermodynamic limit is concentrated in the interval [rn_ , m+] . Since by convexity 
of / we have m + < m(B) for all B > 0, the statements 1) and 2) together imply 
m 2 (0) < m 2 (B) for B > 0. Altogether Eq. (fl3|) is thus established. It is also clear 
that unless the probability distribution is concentrated at the end points of the interval 
[m_,m+] the inequality m 2 (0) < m 2 {B) is strict. This is exactly what happens for 
the Heisenberg ferromagnet because as remarked in |15] the probability distribution is a 
monotonously decreasing function of |m|. 

3. Bosonic Hamiltonians 

We start with the well-known Hamiltonian for bosons in a large box of volume V, 
expressed in terms of the second-quantized creation and annihilation operators a k ,a£ 
satisfying the canonical commutation relations, 

H = fc 2 OkO k + — Y KP) a k+p a q-p a k a q (17) 

k k,p,q 

(with h = 2m = 1). Here, v is the Fourier transform of the two-body potential v(r). We 
assume that there is a bound on the Fourier coefficients |f(k)| < <p < oo. The case of 
hard core potentials can be taken care of by cutting off the potential at some finite value 
that is taken to infinity at the end of the calculations. 

We shall work in the grand canonical ensemble and hence define = H — fiN where 
fi is a chemical potential and N — ^ k a k a k the number operator that generates the 
gauge symmetry. We also add a gauge symmetry breaking term and define 

H^ x = H^ + W(Xa + X*a* ) (18) 

with a complex parameter A that we can take to be a real number without loss of 
generality. While commutes with N this is no longer the case for H^ x if A ^ 0. 
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We now come to the c- number substitution for a and a . Let z be a complex 
number, \z) — exp{— |z| 2 /2 + za^} |0) the coherent state vector in the a Fock space 
and let H(z) — \z)(z\ be the projector onto this vector. There are six relevant operators 
containing a in H^ \, which have the following expectation values [T7] (called lower 
symbols) 

(z\a \z)=z, (z\a a \z) = z 2 , {z\a* a \z) = \z\ 2 

{z\a* \z) = z*, (z\a* a* Q \z) = z* 2 , (z|a a a a |z) = \z\ 4 . 

Each also has an upper symbol, which is a function of z (call it u(z) generically) such 
that an operator F is represented as F — J d z u(z)U(z) , where d 2 z = ir~ 1 dxdy with 
z — x + iy. These symbols are 

a -> z, a a -> z 2 , a* a -> \z\ 2 - 1 

* * * * *2 * * ll4/ill2io 

a — > z , Oo a o ~~ ¥ z > a o a o a o a o ~~ * \ z \ ~~ 4|z| + 2. 

We denote by H' x (z) and H" x {z) the operators obtained from by replacing the 
polynomials in a und (1q by their lower and upper symbols respectively. These operators 
act on the Fock-space of all the modes other than the a mode. One has 

H»{ X )=H'„{z) + S^z) (19) 

with 

W = M + ^ [MM 2 + 2)1/(0) - £ «k<*k(M°) + Kk) + f (-k))] ■ (20) 
The partition functions corresponding to these operators are defined by 



E(fi,X)=Ti H exp[-l3H ll:X } (21) 
e PV P '(^) = s / (/i! x) = J d 2 z Tiw exph/3 ^ x{z)] (22) 

e /3T/ P "(M,A) _ g//^ A) = y d 2 z Trw , e^p^w aW ] (23 ) 

where H is the full Hilbert (Fock) space and Tt' is the Fock space without the a Q mode. 
The functions A), A) and A) are the corresponding finite volume pressures 
and we shall soon sec that in the thermodynamic limit they all coincide. This follows 
from the inequalities 

S'(/i, A) < S( M , A) < E"(/i, A) < E'(n + 2<p/V, X)e^ + ^ (24) 

which we now prove. The first step is 

S( M ,A)>H'(/i,A), (25) 
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which is a consequence of the following two facts: The completeness property of coherent 
states, J d 2 zll(z) — Identity, and 

{z®<j)\e-P H ^\z®(l>) >e^WI«^k«« =e -ma'^*m ) (26) 

where <f> is any normalized vector in TC . This is Jensen's inequality for the expectation 
value of a convex function (like the exponential function) of an operator. To prove 
we take 4> in (j^Hj) to be one of the normalized eigenvectors of H' x (z), in which case 
exp{(0|— (3H' x (z)\(f))} = (<f>\ exp{— 0H' x (z)}\<fr). We then sum over all such eigenvectors 
(for a fixed z) and integrate over z. The left side is then 3(/i, A), while the right side is 

The second inequality is the Berezin-Lieb inequality El [7] 

H(/i ) A)<H"( / x,A). (27) 

Its proof is the following. Let £ Ti denote the complete set of normalized eigenfunc- 
tions of -ff^.A- The partial inner product \'i>j(z)) = (z\§j) is a vector in TC whose square 
norm, given by Cj(z) — j(z)\^ j(z))w , satisfies / d 2 zcj(z) = 1. By using the upper 
symbols, we can write 

(^■1^1^.) = J d?z{t> 3 {z)\H>i. x {z)\^ 3 {z))= J d\{%{z)\H>; hX {z)\%{z)) C] {z), 

where ^'-(z)) is the normalized vector Cj{z)~ 1 / 2 '^! j(z). To compute the trace, we can 
exponentiate this to write S(/Lt, A) as 

^Texpj-/? [ d\c 3 {z){%{z)\Hl x {z)\%{z))\. 

3 

Using Jensen's inequality twice, once for functions and once for expectations as in (j2l)|) , 
S(/x, A) is less than 

J d 2 z Cj (z)exp{(vI/;(z)|-/3^ A (z)|^.(z))} 

i 

<J2j d\c,{z){^{z)\e^{-pH^ x {z)}\^ 3 {z)). 

3 

Since TrU(z) = 1, the last expression can be rewritten as 

J d 2 z 5>,|B(z) ® exp {-PH'; tX (z)} - S"( M , A) 

3 

and (f2"7) is proved. 
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For the last inequality in (|24[) we have to bound 8^(2) in ([20]) . This is easily done in 
terms of the total number operator whose lower symbol is N'(z) = \z\ 2 + X)k^o a k a k- ^ n 
terms of the bound (p on ^(p) 

\5 fl (z)\<2<p(N'(z) + ±)/V + \ f *\. (28) 

Consequently, S"(/x, A) and S'(/z, A) are related by the inequality 

E"(fi, A) < 5?(n + 2<p/V, \y {M+v/v \ (29) 

Closely related to this point is the question of relating H(/i, A) to the maximum value 
of the integrand in (|2"3")1 , which is max z Tin' exp[— f3H' x ( z )] = exp(/?yp max ). This latter 
quantity is often used in discussions of the z substitution problem, e.g., in refs. pi 118] . 
One direction is not hard. It is the inequality (used in ref. [2]) 

E(H, A) > maxTr„, exp[-l3H^ x (z)], (30) 

Z /"J 

and the proof is the same as the proof of ([2!>]) . except that this time we replace the 
completeness relation for the coherent states by the simple inequality Identity > II(z) 
for any fixed number z. 

For the other direction, split the integral defining S"(/i, A) into a part where |z| 2 < £ 
and \z\ 2 > f. Thus, 

E"(fi,X)<(m aS c'Tr w exp{-pH'; jX (z)] + j f d 2 z \z\ 2 Tt-h> expj— A (z)]. (31) 

Dropping the condition \ z\ 2 > £ in the last integral and using \ z\ 2 < N'(z) = N"(z) + 1, wc 
see that the last term on the right side of (I31|) is bounded above by £ -1 S"(/i, X)[Vp"(fj,, A) + 
1], where p"(n, A) denotes the density in the H'' x ensemble. Optimizing over £ leads to 

S"(/i,A)<2[Vy / (/* ) A) + l] maxTr^exp[-/3^' A (z)]. (32) 

Note that p"(p,,X) is order one, since p"(p,X) and A) agree in the thermodynamic 
limit (and are convex in p), and we assume that the density in the original ensemble is 
finite. By (|2SJ|, H'/ iX > H'^ +2(p/v - \p\ - tp/V, and it follows from (Zfl, (32J| and (J3DJ 
that p max agrees with the true pressure p in the thermodynamic limit. Their difference, 
in fact, is at most 0(lnV/V). This is the result obtained by Ginibre in [2] by more 
complicated arguments, under the assumption of super-stability of the interaction, and 
without the explicit error estimates obtained here. 

To summarize the situation so far, we have four expressions for the grand-canonical 
pressure and they are all equal in the thermodynamic limit 

piM, A) = p'(p, A) = p"( M , A) = p max (/i, A) (33) 

when (/j,, A) is not a point at which the density can be infinite. 
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The expectation values (ao a o)^,A an( i ( o)ai,^ are obtained by integrating (|z| 2 — 1) 
and z, respectively, with the weight W fli \(z), given by 

W^ x {z) = E(fi, \)~ lr Tr H > {z\ exp{-/3^, A }|z). 



We will show that for almost every A, the density W^ t \(C,vV ) converges in the ther- 
modynamic limit to a J-function at the point (max = limy^oo z max / VV, where z m ax 
maximizes the partition function Tr^/ exp{— PH^ x( z )}- 
That is, 

^ _1 <«oao>M,A = ^KO^aI 2 = ^kmaxl 2 (34) 

in the thermodynamic limit. This holds for those A where the pressure in the thermo- 
dynamic limit is differentiable; since p(/i, A) is convex (upwards) in A this is true almost 
everywhere. The right and left derivatives exist for every A and hence the quasi average 
lim A ^ + limy^oo ^ _1 |(a ) Mi A| 2 exists. _ 

If we could replace W ll ,\(z) by W fi fi(z)e~ l3X ^( z+z \ the convergence to a (5-function 
would follow from Griffiths' argument in the same way as for the magnetic model in the 
previous section. Because [H, a ] ^ 0, W^.a is not of this product form. However, the 
weight for S"(/i, A), which is 

W';jz) = E"([i, A) _1 Tr-H' exp{-/?fl£ A (z)}, 

does have the right form. In the following we shall show that the two weights are equal 
apart from negligible errors. 

Equality (j33]) holds also for all A, i.e., p(fi, A) = p"(/i, A) = p max (^i, A) in the thermo- 
dynamic limit. In fact, since the upper and lower symbols agree for a and Oq, the error 
estimates above remain unchanged. (Note that since V^Ioq + ag| < S(N +\) + V/5 for 
any 5 > 0, p(/i, A) is finite for all A if it is finite for A = in a small interval around 
/i.) At any point of differentiability with respect to A, Griffiths' Lemma [T5] (see the 
previous section), applied to the partition function S"(/i, A), implies that W[{ X (£y/V) 

converges to a S- function at some point £ on the real axis as V — > oo. (The original 
Griffiths argument can easily be extended to two variables, as we have here. Because 
of radial symmetry, the derivative of the pressure with respect to Im A is zero at any 
non-zero real A.) Moreover, by comparing the derivatives of p" and p max we see that 
£ = limv^oo Zmax/'v/v", since z max /\/V is contained in the interval between the left and 
right derivatives of p max (/i, A) with respect to A. 

We shall now show that the same is true for W^a- To this end, we add another term 
to the Hamiltonian, namely eF = eV J d 2 z H(z)f(zV~ 1 ' 2 ), with e and / real. If f(Q is 
a nice function of two real variables with bounded second derivatives, it is then easy to 
see that the upper and lower symbols of F differ only by a term of order 1 . Namely, for 
some C > independent of zq and V, 



V [d 2 z\(z\z Q )f(f(zV-^)-f(z Q V-^)) 



< C. 
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Hence, in particular, A, e) = p"(fi, A, e) in the TL. Moreover, if f(() — for |£~CI < 6, 
then the pressure is independent of e for |e| small enough (depending only on S). This 
can be seen as follows. We have 

p"(„, A, e) - p"( M , A, 0) = ± l n (e-WV-v*)} (3fi) 

where the last expectation is in the H'^ ensemble at e = 0. The corresponding dis- 
tribution is exponentially localized at z/VV = ( by Griffiths' Lemma, and therefore 
the right side of (|35p goes to zero in the thermodynamic limit for small enough e. 
In particular, the e-derivative of the thermodynamic limit pressure at e = is zero. 
By convexity in e, this implies that the derivative of p at finite volume, given by 
V^ 1 (F) nx = J d 2 z f(zV~ 1 ^ 2 )W Pi ,\(z), goes to zero in the thermodynamic limit. Since / 
was arbitrary, V Jj^_^ >(5 d\ W^x(CVV) — ► as V — > oo. This holds for all 5 > 0, and 
therefore proves the statement. 

Our method also applies to the case when the pressure is not differentiable in A (which 
is the case at A = in the presence of BEC). In this case, the resulting weights W^x and 
W' x need not be S- functions, but as in the argument for |[T3"]) in the previous section 
Griffiths' method implies that they are, for A ^ 0, supported on the real axis between 
the right and left derivative of p and, for A = 0, on a disc (due to the gauge symmetry) 
with radius determined by the right derivative at A = 0. Convexity of the pressure as a 
function of A thus implies that in the thermodynamic limit the supports of the weights 
Wfi^x and W' x for A ^ lie outside of this disc. Hence (rio)\ is monotone increasing in 
A in the thermodynamic limit. In combination with (|34[) this implies in particular that 

Jim i(a5o >^A=o < Hm Jim i(aJ> ) M , A = limlim T?\{a ) n,*\ 2 • (36) 

Hence ((TJ is established. 

We note that by Eq. (f36|) spontaneous symmetry breaking (in the sense that the right 
side of PS)) is not zero) always takes place whenever there is BEC is the usual sense, 
i.e., without explicit gauge breaking (meaning that the left side of (|36p is non-zero). 
Note, however, that a non-vanishing of the right side of (|36|) does not a priori imply a 
non- vanishing of the left side. I.e., it is a priori possible that BEC only shows up after 
introducing an explicit gauge-breaking term to the Hamiltonian. While it is expected on 
physical grounds that positivity of the right side of ([36]) implies positivity of the left side, 
a rigorous proof is lacking, so far. In the example of the Heisenberg magnet discussed 
in the previous section, equality in (|36[) does not generally hold, but still both sides are 
non-vanishing in the same parameter regime. 

To illustrate what could arise mathematically, in principle, consider a weight function 
of the form 

_ ( V 2 -V+l/V for |C| < l/V 

W^ x=0 (W()^wv(0 = { l/V for 1/^ < ICI < 1 (37) 

[ for \z\ >1. 
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This distribution converges for V — > oo to a 5-function at £ = 0, and consequently there 
is no BEC at A = 0. On the other hand, it is easy to see that the weight function 
wv(C) e ~^* V ^ (with an appropriate normalization factor) converges, for any A > 0, to a 
5-function at £ = — 1 as — > oo, and hence there is spontaneous symmetry breaking. 
An open problem for the mathematician is to prove that examples like (|37p do not occur 
in realistic bosonic systems. 
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